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ABSTRACT 



Gradients in the index of refraction caused by temperature fluctuations create optical 
turbulence. This optical turbulence greatly degrades the performance of imaging and laser 
systems operating in the atmosphere. There are several methods of measuring optical 
turbulence in the atmosphere, but few give a high resolution profile of the turbulence at all 
levels displaying the "layers" of turbulence that exist. 

This thesis is an investigation of the use of a high speed temperature probe coupled with 
a VIZ 9000 meteorological system to measure the optical turbulence in the atmosphere from 
a ground base to approximately 25 km. The processing of the data is done by a program 
developed as part of this research. The outputs of the program are r ,, 0 O , and C n : as functions 
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I. INTRODUCTION 



Since the Space Defense Initiative sparked a push to use laser weapons in 
space, many scientists have devoted their attention to the effects of propagating 
a laser beam (high power) through the atmosphere. The problems of 
propagating a laser beam in the atmosphere can be broken into four basic 
categories of detrimental effects [Ref. 1]. They are: absorption and scattering by 
atmospheric aerosols, stimulated Raman scattering, thermal blooming, and 
optical turbulence. It is important to note that all of these effects depend not only 
on the state of the atmosphere but also on the wavelength of the propagating 
beam and the mode of laser operation (CW or pulsed). Significant absorption 
of incident radiation induces thermal blooming when the absorbed energy alters 
the atmospheric index of refraction. At certain wavelengths, such as around 
3.8 [im, atmospheric "windows" exist where absorption is lower. Raman 
scattering is significant only when the peak power is sufficient to induce a 
nonlinear index of refraction in atmospheric nitrogen and oxygen, inducing higher 
order harmonics of the propagating wavelength [Ref. 1], 

At shorter wavelengths atmospheric turbulence and scattering tend to 
dominate the beam propagation problem. Turbulence scrambles the phase of 
a laser beam across the wavefront, thereby spreading the energy over a larger 
solid angle (^radians to milliradians) than would be expected for a diffraction- 
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limited optical system. Scattering from aerosols spreads the energy over very 
large solid angles, on the order of one or more radians, which severely degrades 
collimation and reduces the beam energy that can be delivered to a small target. 

When light propagates through a medium, the local phase of the wave 
depends on the product of the index of refraction of the medium (n) and the path 
length between the source and the point in question. It is known that the index 
of refraction of our atmosphere at a given wavelength is not constant because 
atmospheric density is not constant. Some contributing factors are the decrease 
of density with altitude, differential heating by the sun, turbulent wind flow, and 
convection that causes localized inhomogeneities in the index of refraction, called 
optical turbulence (small variations in n). The index of refraction in air is 

n = 1 + 7 9 xlO" 6 x— (1) 

T 

where P is atmospheric pressure (millibars), T is temperature (kelvin), and the 
humidity term has been ignored since its contribution is so relatively small 
[Ref. 2]. Consequently, light that originates as a plane wave becomes a complex 
distorted wave after propagating a finite distance, as shown in figure 1-1 [Ref. 3], 
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Figure 1-1 : Effect of Atmosphere on 


a Plane Wave 





There are several systems designed to measure optical turbulence. 
However, there is a need for a system that will produce a high resolution, 
accurate vertical profile of the turbulence verses altitude. 

This thesis continues the development of a single point temperature 
thermosonde system to measure the vertical profile of the temperature structure 
parameter. 
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II. BACKGROUND 



A. THEORY 

From the arguments in the previous chapter, it is apparent that analysis of 
a laser system that propagates a beam from ground to space, or vice versa, 
requires an accurate understanding of atmospheric turbulence effects. Although 
these effects are random and complicated, statistical quantities can describe 
many of the relevant phenomena. A successful method that has been developed 
for such analysis uses structure functions to quantify the statistics of atmospheric 
turbulence (rather than autocovariance or autocorrelation functions). Wavefront 
phase distortion depends on localized differential variations in n, rather than upon 
those properties that are common to all portions of the wavefront. Therefore, 
structure parameters are inherently differential in form. A structure tensor as 
defined by Tatarski is [Ref. 2]: 

D ij (T 1 ,T 2 )^<[(f(T 1 )-J{r^)-(f(T 2 )-J{ri))] 2 > (2) 

where the angled brackets represent an ensemble average. 
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By making three assumptions [Ref. 4]: 



1- Isotropy, so that only the magnitude of r is important, not the spatial 
orientation of the measurement 

2- Homogeneity, which implies that a change in position does not alter the 
statistical properties of the turbulence 

3- Incompressible turbulence 

then: 

D f = < [fir,) -f ir 2 ) ] 2 > , ( 3 ) 

and 

r = | ~,-~ 2 | . ( 4 ) 

Furthermore, by the Kolmogrov theory of turbulence, as long as r is within the 
"inertial subrange of turbulence" [Ref. 5], then 

D f = C f 2 r 2/ 3 ( 5 ) 

where C, 2 is a structure parameter. This inertial subrange of turbulence is the 
range from the small scale l 0 of turbulence to the outer scale L 0 . The small scale 
corresponds to the eddy size where viscous dissipation increasingly consumes 
the energy in the eddy. The outer scale corresponds to the largest scale size for 
which eddies can be considered still isotropic (see figure 2-1)[Ref. 4], 
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Figure 2-1: Hierarchy of Eddies for 
Energy Cascade 



For temperature fluctuations that follow a Kolmogrov spectrum the temperature 
structure parameter C T 2 is [Ref. 2] 



2 _ < [r(r 2 ) -rfr,) ] 2 > 

c T - — 



(r 2 -ri) 



2/3 



where T is temperature (kelvin) and r is distance (meters). 

Similarly the index of refraction structure parameter is [Ref. 2] 



(6) 



< [n(r 2 ) -n (r x ) ] 2 > 

' 

where n is the index of refraction and r is distance (meters). Since atmospheric 
turbulence pressure fluctuations are small, the pressure term of the partial 
derivative of equation (1) is negligible, and C n 2 and C T 2 are related by the 
temperature partial derivative [Ref. 6]: 

C 2 = (79xl0* 6 — ) 2 C? ( 8 ) 

rp2 



There are two statistical parameters that characterize the effect of optical 
turbulence. The first is r 0 , which is a coherence length. It is a measure of 
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transverse spatial coherence of wave’s electric field [Ref. 7}. In the absence of 
turbulence, r 0 is comparatively large (see figure 2-2). The importance of r 0 is that 
it describes the spatial resolution of an optical system. It is often called the 
"seeing parameter" or "Fried’s parameter", and it determines the effective 
diffraction limit of an optical system. A large r 0 is desirable if resolution is to be 
determined by the characteristics of the radiating aperture, rather than by 
turbulence. 

In the absence of turbulence, the resolution of a diffraction-limited system 
is 



e c « ^ . (9i 

where 0 D is the angular divergence of the beam (in radians), k is the wavelength. 
D is the aperture size, and a is a geometric factor of order unity. Beam spread 
decreases with increasing aperture size D, until D reaches r 0 . Further increase 
of D beyond r 0 does not result in decreased beam divergence; the system's 
angular resolution stays constant 
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Figure 2-2: Symbolic Projection of r 
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Typical values for r 0 are 5 cm (day) to 10 cm (night) for A. = 500 nm [Ref. 8], For 
a plane wave, r G depends upon the path integral of the index of refraction 
structure parameter, C n 2 . From Ref. 6 



where k is the wave number. Equation 10 shows that r 0 is proportional to X 65 . 

Another related parameter is the isoplanatic angle, 6 0 . Whereas r 0 is the 
spatial coherence length, 0 O is the angular coherence length. It defines the 
"cone" in the atmosphere where turbulence is effectively the same for different 
angular paths. As with r 0 , an optimized optical system should be characterized 
by a large B 0 . 6 0 can be found by two methods. First it can be calculated by 
knowing C n 2 along an optical path [Ref. 9]: 



L 




( 10 ) 
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where z is a path length in meters, or it can be measured by scintillation of 
starlight caused by random variations of local irradiance induced by turbulence 
[Ref. 4). By use of equations 5 and 6, a balloon system with a high speed 
temperature probe can measure a profile of C n 2 from the ground to approximately 
20 km above the earth. This profile can then be used to calculate 9 0 . Comparing 
optically measured 9 0 values to those computed from a balloon system will verify 
the accuracy of these measurements [Ref. 10]. 9 0 also has a a 65 wavelength 
dependence. Typical values are 3-10 ^radians at X = 500 nm [Ref. 8]. 

Both r 0 and 9 0 are measurable parameters that determine the effects of 
atmospheric turbulence that degrade the collimation and spatial coherence of a 
propagating electromagnetic wave. However, r 0 and 9 0 are independent of each 
other. As equation 10 shows, r 0 is equally weighted along the beam path. 9 0 is 
an indicator of higher altitude turbulence, however, as the z 5 ' 3 weighting in 
equation 1 1 suggests. 



B. SYSTEMS REQUIREMENTS 
1. Processing System 

In 1989 the Naval Postgraduate School began development of a single 
point probe thermosonde system to obtain a C T 2 verses altitude profile. The main 
processing component of the system was a VIZ 9000 meteorological processing 
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system that utilizes VIZ 9000 digital Microsonde II thermosondes attached to 
weather balloons. The earlier version sonde (1543-364 Revx6) has 8 channels 
to collect data, but only 4 of these channels are used by the VIZ system for 
meteorological data collection. The remaining unused channels are used by the 
NPS group to measure the temperature of the atmosphere with a high speed 
thermistor probe; however, only 3 of these 4 channels were available for 
measuring temperature due to manufacturer design. Newer version sondes 
(1 543-1 68B) now have all four spare channels available to measure temperature 
with an added probe. 

The earlier NPS single point temperature probe systems used a 0.5 x 
0.5 mm square Thermometries thermistor flake. Since this thermistor proved to 
be too slow, this thesis concentrated on a faster probe [Ref. 10]. A 5 mil VECO 
thermistor probe (page 10) was connected to the thermosonde PC board as 
shown in figure 2-3. 
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Figure 2-3: Connecting Probe To Sonde 
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The VIZ 9000 samples 8 channels every 1.23 ± 0.01 second. These include a 
calibration signal, pressure, temperature, relative humidity and four spare 
channels (1 through 4). With the exception of the pressure channel, all of the 
measurements are resistances. The sampling time of each spare channel varies 
from 32 to 53 ms. The variable time is a consequence of an algorithm built into 
the MICROSONDE II sondes that compensates for the 10 4 to 10 7 range of 
resistance values normally encountered. The resistances measured by the 
system on the unused channels are stored into a computer file at the end of the 
flight and then are converted into temperatures. 

System noise measurements of the older version VIZ microsonde II 
were done by NPS November 1990 [Ref. 10], Table 1 is a summary of the 
measured system noise for the new version sondes. Fixed resistors were 
measured by the sonde giving total system noise expressed in terms of the 
standard deviation of the resistances and equivalent RMS temperatures. 



TABLE 1. VIZ MICROSONDE II SYSTEM NOISE 



Resistor 


o (ch 1) (Q) 


o (ch 2) (Q) 


o (ch 3) (Q) 


a (ch 4) (Q) 


4 $ kQ 


3.9 

=0 005 K 


3.7 

-0.005 K 


3.5 

-0 005 K 


3.6 

-0.005 K 


43 kQ 


12.0 

=0 008 K 


11.8 

-0.008 K 


11.1 

-0 008 K 


10.9 

-0.008 K 


100 kQ 


21.8 

= 0 007 K 


20 7 

-0.007 K 


21.1 

-0 007 K 


22.4 

-0 007 K 


1 M 


345 

-0.012 K 


351 

-0.012 K 


328 

-0.012 K 


333 

-0.012 K 
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43 kQ 


12.0 

=0.008 K 


11.8 

=0.008 K 


11.1 

= 0.008 K 


10.9 

=0.008 K 


100 kQ 


21.8 

=0.007 K 


20.7 

= 0.007 K 


21.1 

= 0.007 K 


22.4 

=0.007 K 


1 M 


345 

=0.012 K 


351 

=0.012 K 


328 

=0.012 K 


333 

=0.012 K 



2. Thermistor Probe 

After trying many different probes, the NPS group has had the most 
success with a 5 mil VECO thermistor probe. Figure 2-4 is a diagram of the 
construction of the probe. The response time of the probe was found to be 



- = A - BjV 



( 12 ) 



where t is the response time (sec), A and B are constants (A = 0.00758 msec' 1 
and B = 0.0077 msec 1 -sec 0 5 /m 05 ), and v is the wind velocity (m/s) across the 
probe. 




3/16" DIA x 12" TUBING 



EPOXY SEALER 

#22 SOLID SILVERED WIRE 

5 mil VECO THERMISTOR 

Figure 2-4: NPS Single Point Probe 
Construction 
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Figure 2-5 shows 1 /t verses v° 5 for wind velocities of 0, 1,2 and 4 m/s measured 
in the lab with a chopped Halogen lamp illuminator as a heat source. The 
response time for the probe for an average ascent rate of 4 m/s was 45 ± 5 
msec. 

To calibrate the probe, a method similar to the Steinhart-Hart equation 

[Ref. 11] 



- = a+b{ In — )+c(ln — ) 2 + d(ln — ) 3 

T r o r o r Q 



( 13 ) 



was found to work best by NPS in June 1990 where a, b, c, and d are 
coefficients and r D room temperature resistance [Ref. 12], The method uses an 
in flight calibration of the probe with the VIZ temperature sensor as a reference. 
A least squares algorithm computes the coefficients that are used to convert the 
thermistor resistances to temperatures. 
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Figure 2-5: Probe Response Time 



The structure function used to calculate C T 2 in equation 6 is sensitive 
to the response time of the probe and the 32 - 53 msec integrating period of the 
A-D converter that is used in the VIZ system. Therefore the structure function 
over a limited scale length becomes [Ref. 11]: 

C(f.T 1 .t 2 )^t/sin» A<r-«[ S ' n( * T ’ /2) fl — '—]dk <”> 

a 2 (kt.,/ 2) 1 +(o>x 2 ) 2 



where a is 2 tt/L 0! the limiting lower frequency (outer scale L 0 ), b is 2 tt/I 0 , the 
limiting upper frequency (inner scale l 0 ), k is 2 tt/A, the wave number, r, is the A-D 
integrating time, is kV where V is the wind velocity (thermosonde ascent rate), 
and t 2 is the time constant of the temperature probe (the upper frequency cutoff). 
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The ratio of the filtered function to the ideal unfiltered structure function 



D(r, 0,0) ' D a 



gives the degradation to the calculation of C T 2 due to high frequency filtration by 
the probe. Figure 2-6 shows the correction required for different separation 
distances between temperature measurements for different probe separations. 



THERMOSONDE PROBE CALIBRATION 
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5 

UJ 

C 02 

0°000 1.000 2.000 3.000 4.000 

PROBE SE PE RATION DISTANCE (m) 

Figure 2-6: Probe Calibration 



3. Computer Program 

The final component of the single point system is a computer program 
written in compiled QuickBASIC 4.5 that uses two converted binary to ASCII VIZ 
flight files as its input. The first file (flight. MTX) provides time, temperature, 
pressure, humidity, spare 1, spare 2, spare 3, and spare 4 resistance data. The 
second file (flight. MDX) provides time and altitude information. 
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A 21 point running mean is applied twice to the altitude data to remove 
the pressure sensor noise. Several different filtering methods were tried, with the 
double 21 point mean giving the best results as shown in figures 2-7 through 
2-11. Next, the resistance data from all four spare channels were converted to 
temperature by the method discussed on page 12. The temperature differences 
from channel 1 to the next data measurement on channel 1 (spare 1-1) were 
calculated along with r (the distance between the measurements) and C n 2 was 
computed according to equation 6. C n 2 was next used to calculate r 0 and G 0 , 
using equations 8 and 9. This same procedure was repeated for temperature 
differences between channels 2 to 1, 3 to 1, and 4 to 1. Summarized, the times 
between temperature measurements where found as follows: 

1) spare 1-1 = T^ - T li+1 (=1.23 sec apart) 

2) spare 2-1 = T 1( - T 2i (=0.05 sec apart) 

3) spare 3-1 = T,, - T 3i (=0.10 sec apart) 

4) spare 4-1 = T tj - T 4i (=0.15 sec apart) 
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Figure 2-7: Effect of Noise in Pressure Sensor & Altitude Data 
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Figure 2-8: Effect of 5 Point Filter of Altitude Data 
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Figure 2-9: Effect of 21 Point Filter of Altitude Data 
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Figure 2-11: Effect of Double 21 Point Filter 
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III. EXPERIMENTS 



A. PROCEDURE 

1. Maui Flights 

The mission of the Maui trip was to aid in site selection for a proposed 
Air Force, 4 meter telescope to be installed. A series of acoustic sounder 
measurements and balloon launches occurred at the Air Force Maui Observation 
Station (AMOS), Haleakala, Hawaii between 26 June to 2 July 1991 [Ref. 13], 
Additionally, the 1.6 meter telescope with the Itek Compensated Imaging System 
(CIS) at AMOS gathered stellar image full width half maximum data. The NPS 
group also launched single point probe thermosonde systems on 28, 29, and 30 
June 1991. The 28 and 29 June flights used the newer (168B) Microsonde II 
sondes, while the 30 June flight used the older version (Rev X-6) sonde. 

All three thermosonde flights were launched after sunset 
(approximately 2200 local time) to avoid solar heating of the probes [Ref. 10]. 
At the time of launch, the NPS group did not know the new version sondes (28 
and 29 June) had different electronic configurations than the older version 
sondes, therefore the probe wiring was connected similarly to the previous 
sondes. This became a problem because a resistance voltage divider and diode 
had been added to the circuit connected across spare channel 4 to monitor the 
sonde battery voltage; the 30 June flight however, being an older version sonde 
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used in previous NPS work, did not include this circuitry. Although the resistance 
divider and diode contaminated the results, particularly at higher altitudes where 
the resistances are in the 100k to 1M range, the lower altitude data was not 
contaminated excessively. 

2. NPS Flights 

The NPS group measured 0 o optically 15 November 1991 from 2130 
to 2330 local time. A thermosonde system was launched from the roof of the 
Naval Postgraduate School’s Spanagel Hall during this same period. The sonde 
was a new version Microsonde II with the voltage divider and diode removed, 
however this launch included two probes as shown in figure 3-1. 
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Figure 3-1: Connecting Probe to New 
Version Sonde 



Two thermistor probes connected to channel 1 and channel 2 respectively and 
separated by 1m provided a fixed separation distance between probes. This 
provided a means to compute horizontal differential temperature measurements. 
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3. MT. Wilson Flights 

In support of characterization measurements of Mt. Wilson, California, 
three balloon launches occurred in September 1991. None of these yielded 
useful data, since RF interference on Mt. Wilson contaminated the measurements. 
The Loran and Omega navigational aid systems used to determine the sondes’ 
location for determination of wind data were particularly suspectable to the 
interference. 

B. RESULTS 

1. Maui Flights 

Table 2 shows the results of the thermosonde flights for 28, 29, and 
30 June 1991 (one flight per night) in comparison with the acoustic sounder data 
and the AMOS 1.6 m telescope data (their measurement of r 0 at the time of the 
balloon flight). 

A 5 point running mean of the temperature data was used in the thermosonde 
C T 2 program and then corrected as shown on page 13. 



TABLE 2. CORRECTED R 0 (cm) COMPARISON MEASUREMENTS 



DATE 


AMOS 


NPS 


1-1 


3-1 


4-1 


AVE 


28 JUN 




10-11 


3.9 


3.1 


3.3 


3.4 


1 29 JUN 


7-8 


9-10 


2.26 


2.25 


2.29 


2.27 


30 JUN 


9-10 




8.9 




7.0 


7.95 
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The flights on the 28 and 29 June differed from the other systems by a factor of 
three. This is believed to be due to the diode circuit in spare channel 4 (new 
version sondes). The 30 June flight (old version sonde) shows agreement with 
the AMOS 1.6 m CIS measurements within 12% discussed in Reference 13. 

A problem exists when computing vertical temperature differences. 
Since the mean atmospheric temperature changes rapidly with altitude vertical 
temperature differences <(T(z 2 ) - T(z,)) 2 > contain some unavoidable mean value 
changes that are not induced by turbulence. This phenomena is particularly 
prominent as the sensors cross a sharp temperature inversion. This thesis 
considered a variety of running means to remove the T 2 - T, mean. Figures 3-2 
and 3-3 show the effect of 5 to 21 point running means on the computation of R 0 
and 0 O . 

Figure 3-2 shows a monotonic increase in r 0 as the length of the 
running mean becomes shorter. This demonstrates that there may be no 
acceptable running mean length and other reduction techniques may be 
necessary. From figure 3-2 it is apparent that the channel 4-1 measurements are 
much less sensitive to the form of the running mean. This is simply a 
consequence of the short time period and distance between measurements. The 
channel 1-1 measurements are widely separated and are inherently sensitive to 
gradients in temperature. In table 2 the channel 4-1 measurements are the most 
accurate. The channel 1-1 measurements in table 2 have been artificially 
enhanced by the 5 point running mean compared to the 21 point results. 
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Another source of systematic errors present when differencing channel 
1 and channels 3 or 4 are the fluctuations in time measurements from channel 
1 to 3 and 1 to 4. To calculate the separation distance r used in equation 10, an 
average of 100 msec and 150 msec time was used in the program. However, the 
signal time occasionally varies by as much as 10% . This leads to noise in the 
separation distance for channels 3-1 and 4-1 that is not present for channel 1-1 
differencing. 

RUNNING MEAN SENSITIVITY 

FLIGHT AF063021 
R_o VS RUNNING MEAN SIZE 

12 

11 . 

LEGEND 
— CHANNEL 1-1 
- CHANNEL 4-1 

7 

6 

5 7 9 11 13 15 17 19 21 

SIZE OF RUNNING MEAN 



Figure 3-2: Running Mean Sensitivity 
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RUNNING MEAN SENSITIVITY 

FLIGHT AF062922 

THETA_o vs running mean size 
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LEGEND 
-CHANNEL 1-1 
-CHANNEL 4-1 



4 

5 7 9 11 13 15 17 19 21 

RUNNING MEAN SIZE 



__ 9 




5 



Figure 3-3: Running Mean Sensitivity 



2. NPS Flight 

Table 3 shows the results of the thermosonde flight done on 15 
November 1 991 compared to optical measurements of the isoplanatic angle. The 
dual point probe and optical measurement method agreed very well, while the 
single point probe are 13% smaller than the optical measurement. 
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TABLE 3. CORRECTED G 0 (^radians) COMPARISON MEASUREMENTS 



DATE 


OPTIC 


DUAL 


1-1 


3-1 


4-1 


AVE 


15 NOV 


3-4 


3.3 


3.1 


2.2 


2.4 


2.6 



Tables 2 and 3 suggest that the sonde circuitry may be inducing 
extraneous noise in the data. The VIZ 9000 MICROSONDE II sondes measure 
resistance by counting clock pulses of an oscillator circuit [Ref. 10]. These clock 
pulses (counts) are then converted to resistance by the manufacturers software 
that is part of the VIZ system. To maintain high resistance and temperature 
resolution (0.002 to 0.008 K), the sonde alters the measurement period by factors 
of two to keep the counts in the 2 15 to 2 16 range. It was found that this altering 
or "ratcheting" happened at approximately the 60 and 260 kQ resistance 
measurements. This has an unfortunate side effect on the small thermistor probe 
self heating. Ohmic heating of the thermistor is proportional to the power 
dissipated by the probe as well as the time duration of the probe current. Figure 
3-4 shows a ratcheting effect in the channel 2 minus channel 1 temperature 
differences at the 60 kQ shift point. A consistent shift of 0.5 K occurs above 7300 
m in the temperature differences as the sonde shifts to a shorter count interval. 
Significant temperature spikes on the order of 0.05 K occur between 7100 and 
7300 m. 
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These fluctuations appear to be the result of an instability near the threshold 
point, as if the sonde shifts periods back and forth several times before 
stabilizing. This effect induces a false turbulence signal in the isoplanatic angle 
on the order of 5-10%. 

This phenomena appears in every sonde launch investigated to date 
when the sonde clock circuitry shifts averaging times. Appendix E includes a 
series of plots that further demonstrates this effect. Surprisingly, the "stutter" 
effect appears in the temperature differences but not in the raw counts recorded 
by the system. It is clear that channels 5-8 each shift periods at different times. 
The cause of this phenomenon is not completely understood at this time and 
should be investigated. 
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FLIGHT AF063021 
ALTITUDE VS DELTA T 
CHANNEL 4-1 




DELTA : (K) 



Figure 3-4: Ratcheting Effect on Delta T 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The new high speed 5 mil VECO thermistor probe proved to be a large 
improvement in the single point probe method of measuring atmospheric 
turbulence. The results of measurements of r 0 and 6 0 are now within 10 to 12 
percent of the other methods of measuring the turbulence. 

The new version sondes are more useful than the older version by allowing 
the use of all 4 spare channels for data collection. The noise in the new version 
sonde is consistent with what was measured in the old version, therefore the 
noise filtering methods used for the old sondes are the same for the new version 
sondes. A successful method of filtering the altitude data was found to be a 
double 21 point running mean. 

This thesis found that a running mean applied to first order temperature 
differences was not sufficient for trend removal. Coherence length computed 
with 5-21 point running means showed a monotonic increase for shorter running 
means. It appears that no running mean is suitable. Other data reduction 
techniques will be necessary. 
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Systematic spikes in Microsonde II temperature data were found near the 
resistance range change points. These temperature spikes contaminated the 
isoplanatic angle by 1 0%. They will have to be removed by software or hardware 
modifications in future work. 

The improved high speed thermistor, along with the new version sonde and 
a computer program incorporating altitude filtering and high order temperature 
differencing will allow the NPS group to have a working single point probe 
system. This system has the ability to plot altitude verses C n 2 displaying where 
the layers of optical turbulence are encountered. 



B. RECOMMENDATIONS 

To gain more confidence in the filtering methods used in the computer 
program and to establish a statistical data base for measurement comparisons, 
a number of flights should be launched with the two point probe system. These 
obviously should cover a wide range of both "good seeing" and "poor seeing" 
days. Another advantage of doing numerous flights to establish a data base is 
that a calibration correction factor could then be applied to the 3-1 and 4-1 
channels. 
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A continued study of how to handle the temperature inversions with a signal 
filter needs to be carried out and incorporated into the computer program. This 
is probably the future research area that would produce the greatest 
improvement to the system. 

It is unknown if the variations in the timing intervals between channel 1 to 
4 can be handled or removed well enough to bring into agreement 
measurements across all channels. If the noise is not improved, then output 
measurements should be restricted to differencing channel 1-1 for a single point 
probe or 2-1 for a dual point probe. 

Shortening the time between scans is a hardware improvement to the VIZ 
Microsonde II that will improve performance. If this scan interval could be 
shortened from the present 1.23 sec to approximately 0.5 sec with the same low 
noise, it would be a big improvement for use as a single point system. This is 
because 0.15 sec tends to be too short for a 50 ms probe response time, while 
1.23 sec is too long (which corresponds to =5 m spacing). Therefore 0.5 sec 
(=2 m spacing) should improve the measurements. 
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APPENDIX A. PROGRAM FOR COMPUTING C n 2 , R 0 , AND Q 0 



DECLARE SUB Crout (a!(), b!(), c!) 

DECLARE FUNCTION temp! (r!, b!()) 

DEFINT I-N 
OPTION BASE 1 
DIM xm ( 0 TO 10) 

DIM a (4 , 4) , b ( 4 ) 

' This program reads the ASCII .mtx and .mdx file created by 
' the VIZ 9000 SYSTEM and performs a least squares fit 
' creating a temperature data output. 

' It then filters the altitude data and computes CN~2, Ro, and 
' theta_o for 4 spare channels of data from the temperature 
' data creating an altitude verses CN^2 ' output file to be 
' plotted with SUPERCALC. 

1 DLW Sep 11 1990 (temperature conversion program) 

' CRH Nov 10 1991 

1 To declare dynamic data array, invoke Quickbasic with /ah 
' option 



' LIST 
' start 
' stp = 

' dta ( i 
' dta ( i 



OF VARIABLES AND CONSTANTS: 

= starting altitude data point 
stopping altitude data point 



dta ( i 
dta ( i 

dta ( i 
dta ( i 
dta ( i 
dta ( i 
dta ( i 
dta ( i 
dta ( i 
dta ( i 
dta ( i 
dta ( i 
dta ( i 
dta ( i 
dta ( i 
b ( 1 ) = 

b ( 2 ) = 



, 1 ) = 

/ 2 ) = 

, 3) = 

, 4 ) = 



time (sec) 
pressure (mbar) 



VIZ temp 
humidity 



(C) 

(%) 



,5) = 
, 6 ) = 
,7) = 
, 8 ) = 
, 9 ) = 
, 10 ) = 
,H) = 
, 12 ) = 
,13) = 
,14) = 
,15) = 
,16) = 
, 17 ) = 



spare 

spare 

spare 

spare 



1 

2 

3 

4 



temp 

temp 

temp 

temp 



rol= Ro for spare 1-1 
Ool= theta_o for 
spare 1-1 

ro2= Ro for spare 2-1 
Oo2= theta_o for 
spare 2-1 

ro3= Ro for spare 3-1 
Oo3= theta_o for 
spare 3-1 
(K) 

(K) 

(K) 

(K) 



resistance 
resistance 
resistance 
resistance 
VIZ altitude (m) 
first 21 pt mean altitude 
2nd 21 pt mean altitude 
delta T for all spares (reused) 

delta T for all spares (reused) 
for spare 1-1 ro4= Ro for spare 

for spare 2-1 Oo4= theta_o for 

for spare 3-1 spare 4-1 

_ for spare 4-1 
coefficient of polynomial to convert resistance to 
temperature 



mean 
CN~2 
CN A 2 
CN A 2 
CN~2 



4-1 
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' b(3)= coefficient of polynomial to convert resistance to 
temperature 
' b(4)= , , 

' low = number of data pts lost due to running mean size 
' high = number of data pts lost due to running mean size at 
end of flight for spare 1-1 only. 

' ST = desired starting altitude 
' RM = desired running mean size 
' pts = number of vertical data points 



REM $DYNAMIC 
DIM dta ( 6000 , 17) 

tzero = 273.15 ' zero Celcius 

INPUT "enter MTX filename MTX$ 

PRINT 

PRINT MTX$ 

OPEN MTX$ FOR INPUT AS #1 

INPUT "enter MDX filename MDX$ 

PRINT 

PRINT MDX$ 

OPEN MDX$ FOR INPUT AS #2 
1 = 0 

DO UNTIL EOF ( 1 ) 

1 = 1 + 1 

INPUT #1, dta ( 1 , 1), dta ( 1 , 2), dta(l, 3), dta ( 1 , 4), 

dta ( 1 , 5), dta(l, 6), dta ( 1 , 7), dta ( 1 , 8) 

LOOP 

1 = 0 

DO UNTIL EOF (2) 

1 = 1 + 1 

INPUT #2, x, y, dta ( 1 , 9) 'x and y are unneeded variables 
LOOP 

•'FIND START AND STOP OF FLIGHT 
N = 1 



FOR i = 


1 


TO 


(N 


/ 


2) 








IF (dta ( i , 


9 ) 


< 


dta ( i 


+ 1 / 


9 ) ) 


AND 


(dta ( i 


+ 


i, 


9 ) 


< 


dta ( i 


+ 2, 


9 ) ) 


AND 


( dta ( i 


+ 


2, 


9 ) 


< 


dta ( i 


+ 3, 


9 ) ) 


AND 


(dta ( i 


+ 


3, 


9 ) 


< 


dta ( i 


+ 4, 


9 ) ) 


AND 


(dta ( i 


+ 


1, 


9 ) 




dta ( i , 


9 )) 


> 2 


AND 
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dta ( i +1, 9)) >2 



(dta ( i + 2, 9) - 

THEN GOTO 10 
NEXT i 

10 start = i 'exited loop when conditions were met 

FOR i = N TO (N / 2) STEP -1 

IF dta ( i , 9) > dta ( i - 1, 9) AND 

dta ( i - 1, 9) > dta ( i - 2, 9) AND 
dta ( i -2,9) > dta ( i - 3, 9) THEN GOTO 20 
NEXT i 

20 stp = i 'exited loop when conditions were met 

CLOSE #1: CLOSE #2 

'CORRECT ALTITUDE (VIZ MAX ALTITUDE INPUT IS 2500 m) 

INPUT "enter altitude correction (m)"; zcor 
PRINT 

PRINT zcor 

FOR i = start TO stp 

dta(i, 9) = dta(i, 9) + zcor 
NEXT i 

'DO DOUBLE 21 PT RUNNING MEAN ON ALTITUDE (Z) 

'first 21 point mean 
FOR i = start TO (stp - 20) 
num ! = 0 ! 

FOR j=iTO(i+20) 

num! = num! + dta ( j , 9) 

NEXT j 

dta ( i + 10, 10) = num! / 21! 

NEXT i 

'second 21 point mean 
FOR i = (start + 10) TO (stp - 30) 
num ! = 0 ! 

FOR j = i TO (i+20) 

num! = num! + dta(j, 10) 

NEXT j 

dta(i + 10, 11) = num! / 21! 

NEXT i 

' compute average ascent distance between data points at start 
' of flight 

zlave! = (dta(start + 24, 11) - dta(start + 20, 11)) / 4! 
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' compute average ascent distance between data points at end 
' of flight 

z2ave! = (dta(stp - 20, 11) - dta ( stp - 24, 11)) / 4! 

' compute averaged altitudes to fill in first 20 data pts lost 
' by double filter 

FOR j = 1 TO 20 

dta(j, 11) = dta(start + 20, 11) - ((21 - j) * (zlave!) ) 

NEXT j 

j = 1 

FOR i = start TO (start + 19) 
dta(i, 11) = dta(j, 11) 
j = j + 1 
NEXT i 

' compute averaged altitudes to fill in last 20 data pts lost 
' by double filter 

FOR j = 1 TO 20 

dta ( j , 11) = dta ( stp - 20, 11) + (j * z2ave!) 

NEXT j 

j = 1 

FOR i = (stp - 19) TO stp 
dta ( i , 11) = dta ( j , 11) 

j = j -• 1 

NEXT i 

' The following is a program to convert resistance data to 
' temp data 

' begin least squares fit 
m = 4 

FOR i = 1 TO m 

FOR j = 1 TO m 

a ( i , j) = 0! 

NEXT j 
b ( i ) = 0! 

NEXT i 

m2 = m + m - 2 
FOR k = start TO stp 
xk = LOG (dta (k , 5) ) 

yk = 1 ! / (dta(k, 3) + tzero) 

xm ( 0 ) = 1 ! 
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' xm 



x~m 



FOR i = 1 TO m2 

xm ( i ) = xk * xm ( i - 1 ) 

NEXT i 

FOR i = 1 TO m 
i2 = i - 2 
FOR j = i TO m 

a(i, j) = a(i, j) + xm(i2 + j) ' x A (i-l+j-1) 

NEXT j 

b(i) = b(i) + yk * xm(i - 1) 

NEXT i 

NEXT k 

' set lower half of matrix 

FOR i = 1 TO m 
FOR j = 1 TO i 

a(i, j) = a ( j , i) 

NEXT j 
NEXT i 

c = m 

CALL Crout(a(), b(), c) 

PRINT b ( 1 ) , b ( 2 ) , b ( 3 ) , b(4) 

' convert to temperature 

FOR i = start TO stp 
dta(i, 3) = dta(i, 3) + tzero 
resistance = dta ( i , 5) 
dta ( i , 5) = temp (resistance, b()) 

resistance = dta ( i , 6) 
dta(i, 6) = temp (resistance, b()) 
resistance = dta ( i , 7) 
dta ( i , 7) = temp (resistance, b()) 
resistance = dta ( i , 8) 
dta ( i , 8) = temp (resistance , b()) 

NEXT i 

'RENUMBER ARRAY FROM 1 TO N 

FOR i = start TO stp 
FOR j = 1 TO 11 

dta(i + 1 - start, j) = dta ( i , j) 

NEXT j 
NEXT i 

N = (i - start + 1) 

PRINT 
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PRINT "starting altitude = "; dta(l, 11) 

PRINT "stopping altitude = "; dta(N, 11) 

' allow user to select new starting altitude to examine 
1 portion of flight 

PRINT 

INPUT "enter new desired starting altitude"; ST 

PRINT 

PRINT ST 

FOR i = 1 TO N 

IF dta ( i , 11) > ST THEN GOTO 30 

NEXT i 

30 start = i 
PRINT 

PRINT "New starting altitude = "; dta (start, 11) 

' compute number of data points 
stp = N 

pts = stp - start + 1 
PRINT 

PRINT "number of data points are "; pts 
'ENTER RUNNING MEAN TO BE USED FOR dT 

INPUT "enter size of running mean (odd number) for temp"; RM 
PRINT 

PRINT "dT running mean is "; RM 

'COMPUTE CT"2 FOR SPARE 1 - SPARE 1 TEMP DIFFERENCE 

'COMPUTE DELTA T 

FOR i = start TO (stp - 1) 

dta ( i , 12) = dta ( i , 5) - dta(i + 1, 5) 

NEXT i 

'DO RUNNING MEAN OF dT 

low = (RM - 1) / 2 

high = (RM + 1) / 2 

FOR i = start TO (stp - RM) 
num! = 0! 

FOR j = i TO (i + RM - 1) 

num! = num! + dta(j, 12) 

NEXT j 

dta ( i + low, 13) = num! / RM! 

NEXT i 
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'compute cn2 

FOR i = (start + low) TO (stp - high) 

dta ( i , 14) = (((.000079 * dta ( i , 2)) / (dta(i, 3) A 2)) A 2)* 

( (dta ( i , 12) - dta ( i , 13)) A 2) / ((dta(i + 1, 11) 

- dta ( i , 11) ) A ( . 6667) ) 

NEXT i 

h! = 2.1 * (1.46 * (2 * 3.14159 / 500E-09) A 2) A -0.6 

sum! = 0! 

' compute Ro (CN A 2 * dz) 

FOR i = (start + low) TO (stp - high) 

sum! = sum! + ((dta(i + 1, 11) - dta ( i , 11)) * dta ( i , 14)) 

' dz=dta ( i+1 , 11) -dta ( i , 11) 

NEXT i 

sum2 ! = sum! A (-.6) 
rol! = h! * sum2 ! 

'COMPUTE Oo for 1-1 (CN A 2 * dz * z A 5/3) 

f! = (2.905 * (2 * 3.14159 / 500E-09) A 2) A -0.6 

sumO! = 0! 

FOR i = (start + low) TO (stp - high) 

sumO! = sumO! + (((dta(i +1, 11) - dta ( i , 11)) * 
dta(i, 14)) * (dta(i, 11) - dta(start, 11)) A 1.6667) 

NEXT i 

sum02 ! = sumO! A (-.6) 

Ool! = f! * sum02 ! 

'COMPUTE CT A 2 FOR SPARE 2 - SPARE 1 TEMP DIFFERENCE 

'compute delta T 

FOR i = start TO stp 

dta(i, 12) = dta(i, 5) - dta(i, 6) 

NEXT i 

'do running mean of dT 

FOR i = start TO (stp - RM + 1) 
num! = 0! 

FOR j=iTO(i+RM-l) 

num! = num! + dta(j, 12) 

NEXT j 

dta(i + low, 13) = num! / RM! 

NEXT i 

'compute cn2 
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FOR i = (start + low) TO (stp - low) 

dta ( i , 15) = (((.000079 * dta ( i , 2)) / (dta(i, 3) * 2)) "2)* 
( ( dta ( i , 12) - dta ( i , 13)) ~ 2) / ((((dta(i + 1, 11) 

- dta ( i , 11)) / 1.23) * .05) ~ (.6667)) 

NEXT i 
sum! - 0! 
sum2 ! = 0! 

FOR i = (start + low) TO (stp - low) 

sum! = sum! + (dta(i, 15) * (dta(i + 1, 11) - dta(i, 11))) 

NEXT i 

sum2 ! = sum! ~ (-.6) 
ro2! = h! * sum2 ! 

'compute Oo for 2-1 

sumO! = 0! 

FOR i = (start + low) TO (stp - low) 

sumO! = sumO ! + (((dta(i + 1, 11) - dta(i, 11)) * 

dta ( i , 15)) * (dta(i, 11) - dta(start, 11)) ~ 1.6667) 

NEXT i 

sum02 ! = sumO! A (-.6) 

Oo2 ! = f ! * sum02 ! 

'COMPUTE CT~ 2 FOR SPARE 3 - SPARE 1 TEMP DIFFERENCE 

' compute dT 

FOR i = start TO stp 

dta(i, 12) = dta(i, 5) - dta ( i , 7) 

NEXT i 



'compute running mean of dT 

FOR i = start TO (stp - RM + 1) 
num! = 0! 

FOR j = i TO (i + RM - 1) 

num! = num! + dta ( j , 12) 

NEXT j 

dta(i + low, 13) = num! / RM! 

NEXT i 

'compute cn2 

FOR i = (start + low) TO (stp - low) 

dta ( i , 16) = (((.000079 * dta ( i , 2)) / (dta(i, 3) * 2)) "2)* 

( (dta ( i , 12) - dta ( i , 13)) ~ 2) / ((((dta(i + 1, 11) 

- dta ( i , 11)) / 1.23) * .1) ~ (.6667)) 'dz/1.23= vel'y 

'r=vel'y * O.lsec 

NEXT i 
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sum! = 0! 
sum2 ! = 0! 

FOR i = (start + low) TO (stp - low) 
sum! = sum! + (dta(i, 16) * (dta(i + 1, 11) - dta ( i , 11))) 
NEXT i 

sum2 ! = sum! A (-.6) 
ro3 ! = h! * sum2 ! 

'compute Oo for 3-1 

sumO! = 0! 

FOR i = (start + low) TO (stp - low) 

sumO! = sumO! + (((dta(i + 1, 11) - dta(i, 11)) * 

dta(i, 16)) * (dta(i, 11) - dta(start, 11)) A 1.6667) 

NEXT i 

sum02 ! = sumO! * (-.6) 

Oo3 ! = f ! * sum02 ! 

'COMPUTE CT A 2 FOR SPARE 4 - SPARE 1 TEMP DIFFERENCE 

'compute delta T 

FOR i = start TO stp 

dta(i, 12) = dta(i, 5) - dta(i, 8) 

NEXT i 

FOR i = start TO (stp - RM + 1) 
num! = 0! 

FOR j = i TO (i+RM-1) 

num! = num! + dta(j, 12) 

NEXT j 

dta ( i + low, 13) = num! / RM! 

NEXT i 

' compute cn2 

FOR i = (start + low) TO (stp - low) 

dta ( i , 17) = (((.000079 * dta(i, 2)) / (dta(i, 3) A 2)) A 2)* 
( (dta ( i , 12) - dta ( i , 13)) - 2) / ((((dta(i + 1, 11) 

- dta ( i , 11)) / 1.23) * .15) " (.6667)) 

NEXT i 
sum! = 0! 
sum2 ! = 0! 

FOR i = (start + low) TO (stp - low) 
sum! = sum! + dta ( i , 17) * (dta(i + 1, 11) - dta(i, 11)) 

NEXT i 

sum2 ! = sum! A (-.6) 
ro4 ! = h! * sum2 ! 

'COMPUTE Oo for 4-1 
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sumO! = 01 

FOR i = (start + low) TO (stp - low) 

sumOl = sumO! + (((dta(i +1, 11) - dta(i, 11)) * 

dta ( i , 17)) * (dta ( i , 11) - dta(start, 11)) ~ 1.6667) 

NEXT i 

sum02 ! = sumO! A (-.6) 

Oo4 ! = f ! * sum02 ! 

'OUTPUT TO A FILE ALTITUDE VERSES CN2 FOR EACH CHANNEL TO BE 
'PLOTTED WITH SUPERCALC 

INPUT "enter output file name ", output$ 

OPEN output? FOR OUTPUT AS #3 

FOR i = start TO stp 

WRITE #3, dta ( i , 11), dta ( i , 14), dta(i, 15), dta ( i , 16), 

dta ( i , 17 ) 

NEXT i 



PRINT 



PRINT 


"Roll= " 


, rol ! 


PRINT 


"Ro2= ", 


ro2 ! 


PRINT 


"Ro3= ", 


ro3 ! 


PRINT 


" Ro4= ", 


ro4 ! 


PRINT 


"Ool= " ; 


Ool ! 


PRINT 


"Oo2= " ; 


Oo2 ! 


PRINT 


"Oo3= " ; 


Oo3 ! 


PRINT 


"Oo4= 


Oo4 ! 



END 



REM $STATIC 
DEFSNG I-N 

SUB Crout (a ( ) , b() , c) 

DEFINT I -M 

' this subroutine solves the set of simultaneous equations 
' A*X=B where A is a symmetric matrix using the Crout or 
' Cholesky procedure. 

' The original A and B arrays are destroyed, 
m = c 

all = a ( 1 , 1) 

IF all = 0 THEN GOTO singular 
FOR i = 2 TO m 

a ( 1 , i) = a ( i , 1) / all 

NEXT i 

b ( 1 ) = b ( 1 ) / all 

I 
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a(k, j) 



'set up auxiliary matrix 

I 

FOR j = 2 TO m 
jl = j - 1 
FOR i = j TO m 
asave = 0! 

FOR k = 1 TO jl 

asave = asave + a(i, k) * 

NEXT k 

a(i, j) =a(i ; j) - asave 
IF i > j THEN a ( j , i) = a(i, j) / a(j, j) 
NEXT i 
bsav = 0! 

FOR k = 1 TO jl 

bsav = bsav + a(j, k) * b(k) 

NEXT k 

ajj = a(j, j) 

IF ajj = 0 THEN GOTO singular 
b(j) = (b(j) - bsav) / ajj 
NEXT j 

I 

' back substitution and calculate determinant 

i 

ml = m - 1 
FOR i = 1 TO ml 
bsav - 0! 
mi = m - i 
mil = mi + 1 
FOR j = mil TO m 

bsav = bsav + a(mi, j) * b(j) 

NEXT j 

b(mi) = b(mi) - bsav 
NEXT i 

GOTO finished 
singular : 

PRINT "***Matrix is singular***" 
finished : 

END SUB 

DEFSNG I-M 

FUNCTION temp (r, b()) 
lnr = LOG (r ) 

temp = l! / (b(l) + lnr * (b(2) + lnr * (b(3) + b(4 

END FUNCTION 



* lnr) ) ) 
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APPENDIX B. NOTES ON OPERATION OF VIZ-9000 SYSTEM 



The following are notes to support the operation of the VIZ system. They 
are to be used in conjunction with the manufacturer’s manual. These notes are 
a brief summary of the steps the NPS group found to be important upon initial 
indoctrination of how to use the VIZ system. 

PART I 



• The initial wiring of the system should be as shown in figure B-1. 



COMPUTER 




Figure B-1: Wiring Diagram 
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PART II 



• After properly connecting all the cables, and properly installing the 

navigation and receiver antennas, the system should be booted up ensuring the 

"rack" is turned on before the computer. 

(WARNING: be sure the nav antenna is not too close to the 
receiver antenna before starting the system. The 
receiver antenna makes a complete rotation upon 
energizing and can damage a nearby nav antenna.) 

• The next step is to run "SETUP". SETUP consists of 8 parts. The 

following configuration should be used for the Monterey Bay area: 

★★★Whenever data is selected or entered, to leave the present screen 
and return to the previous menu, press "escape" and follow the 
directions on the screen. 

1) Station Installation: 

Customer Name: NPS 

Station ID: NPS 

LAT: 36° 35’ 42.0 N 

LONG: 121° 52’ 57.0 W 

WMO station no: 0 

Station Ht: 38 m (2500 m MAX) 

Local time: 0 
Surface data timer: 0 
Sonde type: 2 
PTU selection: Y 

2) Receiver Installation: 

start-stop: 395-410 hz 

3) Navaid Installation: 

Zeemet system: LORAN 
Tracker 1: LORAN 
Tracker 2: OMEGA 

****lf desire an omega launch, then just toggle the zeemet system 
to omega (by using space bar), nothing else has to be 
changed. 
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****For cross chain loran flight, change tracker 2 to loran. 



After entering "LORAN, LORAN, OMEGA", go to tracker 1 and press "enter". This 
will open a sub-window to choose what loran chain you want. Highlight "9940" 
and press the space bar so the "check mark" moves to 9940. Next press 
"escape", go to tracker 2, press "enter", and now select for omega: 



North Dakota: common 

La Reunion: include 

Argentina: include 

Australia: include 

4) Wind Computation (done in height) 



2. 1000 90 3 15 

3. 10000 120 4 15 

OMEGA 

1.0 120 4 15 

2. 1000 180 5 15 

3. 10000 240 7 15 

5) Loran Chain File: 

-leave untouched 

6) Antenna Configuration: 



The cycle time needs to be rechecked at least yearly. Presently the cycle time 
is 48.3 sec. 

****Remember at the end of this menu to adjust antenna center for each 
site (need magnetic variation for site). For Monterey, magnetic 
variation is 16° East, therefore compass heading should read 074 
(so yagi antenna points 090 T). 



OMEGA 



Japan: 

Norway: 

Liberia: 

Hawaii: 



include 

exclude 

exclude 

include 



LORAN 

1.0 60 2 15 
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6) Antenna Configuration (continued): 
mode: auto 
select: omm 
systems: both 

switch: 4000 m down range 

When this information is entered, antenna will cycle and stop at its midpoint. 
This is when the orientation should be adjusted (adjust antenna to most probable 
wind azimuth). 

PART III 



• To run a flight: 

1) choose which type of sonde (loran or omega) 

2) fill balloon and attach sonde 

3) get permission from Monterey Airport tower to launch (375-1214) 

4) type "START" 

• When menu appears, enter sonde serial number and flight ID: 

flight ID = month + day + hour (2 digits each) 

EX: 031121 = March 11, 2100 

• Once this is entered, check Navaid status. Use instruments to enter surface 
data. 

• Once sonde is prepared, check receiver status (signal strength should be 200- 
255). Receiver should lock on sonde after a few minutes. If not, go to receiver 
display [Alt-s] to observe sweep of spectrum. 

• Next check calibration. Check to see that temperature and humidity lock- 
in numbers match those on sonde sticker (if close, go with computer, if way off, 
go with sonde label). Enter offsets if needed. 
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• Finally, when "?'s" on 1-4 have changed to arm for launch. You can 
launch the sonde at anytime after you press arm for launch. However, do not 
hold the sonde too long because you only have two hours of allotted time from 
the "arm for launch" time for a complete flight. 

• Wait a few minutes for PTU data to show up on the screen. 

• Once you have PTU data, go to "Data Display Mode". Choose the PTU 
plot display. 

• Check ascent rate for 180-240 m/min. 

• At any time in the flight you can check on the pointing of the antenna and 
S/N ratio levels by going to "utilities". 

****Do not go to receiver mode in utilities once launch has started. This 
will cause a stopping of data collection. 

• When flight has finished, you should be in either PTU plot or general flight 
menu screen. Now go to utilities to save flight data to disk . 

• After this is done, you can exit flight menu. 
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APPENDIX C. PLOTS OF C n 2 



FLIGHT AF062822 
ALTITUDE VS CV2 
Ro = 19, TIIETA_o = 6.7 

CHANNEL 1-1, 5 FT MEAN 




2.2 c 'e-14 



3*-u 



Cn*2 



Figure C-1: Altitude verses C n 2 
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FLIGHT AF0GR!L , ‘2 
ATTITI'DK VS (V2 
Ro = R.;i, THETA 0 = ii.:i 

CHANNEL 1-1, 5 FT MEAN 






Figure C-2: Altitude verses C n 2 
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FLIGHT AF0G30H1 
ALTITUDE VS Cn-2 
Ho = 8.9, THETA.o = 5.4 

CHANNEL 1-1, 5 FT MEAN 



: 



15000 






5M5 If— 14 1.5e-M 2<r— 14 

Crt'2 



Figure C-3: Altitude verses C n 2 
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FLIGHT NSlII.VTi 
ALTITCDK VS (Tr2 






Ho = 1.9, THETA o = 3.3 






CHANNEL 2-1, TWO POINT PROBE 




u 






I 




£ 


F 




1 

U I 






§ : 


h- 




1- 1 


k 




5 1 


i 




< 






















Cn-2 



Figure C-4: Altitude verses C n 2 
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